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Abstract 

The linear complementarity problem, LCP{q, M), is defined as follows. For given M G 
■j^rnxm^q G 7^'^, find z such that q + Mz > 0, z > 0, z^{q + Mz) = 0, or certify that 
there is no such z. It is well known that the problem of finding a Nash equilibrium for a bimatrix 
game (2-NASH) can be formulated as a linear complementarity problem (LCP). In addition, 2- 
NASH is known to be complete in the complexity class PPAD (Polynomial-time Parity Argument 
Directed). However, the ingeniously constructed reduction (which is designed for any PPAD prob- 
lem) is very complicated, so while of great theoretical significance, it is not practical for actually 
solving an LCP via 2-NASH, and it may not provide the potential insight that can be gained from 
studying the game obtained from a problem formulated as an LCP (e.g. market equilibrium). 
The main goal of this paper is the construction of a simple explicit reduction of any LCP{q, M) 
that can be verified as belonging to WAV via the graph induced by the generic Lemke algorithm 
with some positive covering vector d, to a symmetric 2-NASH. In particular, any endpoint of this 
graph (with the exception of the initial point of the algorithm) corresponds to either a solution 
or to a so-called secondary ray. Thus, an LCP problem is verified as belonging to WAD if any 
secondary ray can be used to construct, in polynomial time, a certificate that there is no solution 
to the problem. We achieve our goal by showing that for any M, q and positive d satisfying a 
certain nondegeneracy assumption with respect to M, we can simply and directly construct a 
symmetric 2-NASH whose Nash equilibria correspond one-to-one to the end points of the graph 
induced by LCP{q, M) and the Lemke algorithm with a covering vector d. We note that for a 
given M the reduction works for all positive d with the exception of a subset of measure 0. 



1 Introduction 



The linear complementarity problem LCP{q, M) is defined as 

For given q G n"^, M G 7?."'X"^, find z G TZ"^ such that q + Mz>0, z>0, z^{q + Mz) = 0. 

The LCP is notable for its wide range of applications, from well understood and relatively easy to 
solve problems, such as linear and convex quadratic programming problems, to MV-hard problems. A 
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major effort in LCP theory fiad been tlie study of variants of Lemke's algorithm, a Simplex-hke vertex 
following algorithm. In particular, for a given positive covering vector d, the Lemke((i) algorithm 
goes through a path of adjacent vertices of the 'extended' LCP{q, M) (denoted by ELCP{d,q, M)) 
where d is attached to M with an artificial variable zq. Assuming (without loss of generality) that 
ELCP{d,q, M) is nondegenerate, Lemke{d) is guaranteed to terminate in a finite number of steps 
with either a solution to the original problem or with a secondary ray of ELCP{d,q, M). If the 
secondary rays can certify (in polynomial time) that there is no solution to LCP{q, M), we say that 
the problem is Lemke(d)-resolvahle. One of the major themes of LCP research over the years has 
been the search for classes of matrices M and covering vectors d for which LCP{q, M) is Lemke((i)- 
resolvable for all q. Several such classes (usually applicable for all d > 0) were identified (see e.g. 
|CPS92] . |Mur88] and the references therein). 

The introduction of the WAD [Polynomial-time Parity Argument Directed) complexity class 
in |Pap94| provides an effective and elegant framework for analyzing the complexity of Lemke((i)- 
resolvable linear complementarity problems since, in general, the directed graph induced by the 
Lemke(d) algorithm for a given LCP{q, M) can be used to verify the membership of the problem 
in WAD. We say in this case that the problem is Lemke(d) W AD -verified. This development is 
significant with respect to LCP theory since it has been shown in [ MP91] that if WAD is TVP-hard 
then MV = CoMV, lending support to the long standing informal belief that LCPs resolvable by 
Lemke(c?) algorithm are in some way special. 

What makes the class WAD particularly interesting is the fact that several well known problems, 
such as finding a Brouwer fixed-point, were identified in |Pap94| as T^T-'^P-complete. The discovery, 
in a string of papers ([DP05]. [DP05a| . |CD05] and [CDOSaj ). that finding a Nash equilibrium of 
a bimatrix game (2-NASH) is T-'T^^X'-complete has significant consequences in the context of LCP 
theory. It has been known since the early days of LCP research that the 2-NASH problem can 
be formulated as an LCP with roughly the same size and with the coefficient matrix belonging 
to one of several well known classes resolvable by Lemke((i) algorithm. The fact that 2-NASH is 
WAD-complete means that any LCP{q, M) verifiable as a member in WAD (including all classes 
that contain 2-NASH) can be reduced to a 2-NASH problem. However, the known reduction is 
quite complicated. It requires several stages that involve reducing the given LCP{q, M) to finding 
an approximate Brouwer fixed point of an appropriate function, followed by reducing the latter to 
3-graphical NASH (using small polymatrix games to simulate the computation of certain simple 
arithmetic operations), and finally, reducing the 3-graphical NASH to 2-NASlf]- While there seems 
to be no discussion in the vast literature on LCP suggesting the possibility that Lemke(d) WAD- 
verified LCPs can be reduced to 2-NASH, the discovery that 2-NASH is "PP^lP-complete motivated 
us to search for the existence of a direct simple reduction of such problems to 2-NASH. 

The main result of this paper is the introduction of a direct, simple reduction of almost any 
Lemke((i) T^'P^X'-verified linear complementarity problems to a symmetric 2-NASH. In fact, we 
introduce a stronger result as follows. Consider a generic Lemke((i) LCP{q, M) (which we call 
LLCP{d, q, M)) whose 'solutions' are defined to be either actual solutions of LCP{q, M) or secondary 
rays of ELCP{d,q, M). Obviously this probleirH is Lemke((i) WAD-Yevi&ed. Through a series of 
steps we show how to construct a symmetric bimatrix game whose equilibria correspond one-to-one 

clear 'bird's-eye view' description of the reduction can be found in [DGP09| . 
^ Where, as we assume without loss of generality, its extended form, ELCP{d, q, M), is nondegenerate. 
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to the 'solutions' of LLCP{d,q,M). The point is that if LCP{q,M) is Lemke{d) VVAV-venG.ed, 
the Nash equihbria of the constructed bimatrix game correspond to solutions (or certifications for 
infeasibility) of the given problem, which means that the constructed bimatrix game properly resolves 
an LCP{q, M) if it is Lemke(d)-verified. 

We begin by reviewing the Lemke((i) algorithm (in Section [2]) and bimatrix games (in Section 
[3]). Next, we introduce (in Section H]) the complexity class WAT), and briefly discuss the no- 
tion of Lemke((i) WAV-venHed LCPs. In addition, we present the majority of matrix classes 
known to be Lemke((i) VPA'D-veiihed, and identify a number of matrix classes whose corresponding 
LCPs are WAV-complete. We conclude the introductory sections by introducing (in Section [5]) 
LLCP{d,q,M), the generic Lemke{d) LCP{q,M). 

Our main results are presented in sections [GHSl We start by introducing in Section [6] a very simple 
reduction of LCP{q, M), where M belongs to a class of matrices for which a solution is guaranteed to 
exist for all q, to a symmetric 2-NASH. The cost matrix of the resulting bimatrix game is composed of 
M with an extra row and column. In particular, we show that the solutions of the given LCP{q, M) 
correspond one-to-one to the Nash equilibria which use with positive probability for the pure strategy 
corresponding to the extra column of the cost matrix. Moreover, we show that the Nash equilibria 
which do not use the pure strategy corresponding to the extra column of the cost matrix of the 
resulting game correspond one-to-one to the so-called 'secondary directions' of ELCP{e, q, M). Note 
that at this stage we address only e - the vector of all ones - as a covering vector, and that we do not 
reach yet our goal as the reduction may produce secondary directions rather than secondary raysH 

In Section [71 we extend the basic reduction above (by considering an augmented problem) so 
that the constructed bimatrix game produces either a solution for LCP{q, M), a secondary ray for 
ELCP{e, q, M), or a non-zero vector which is a solution to LCP{e, M) (and is actually also a special 
case of a secondary direction of ELCP{e, q, M)). 

In Section [H we show that if a secondary direction generated by the bimatrix game constructed 
in the previous section is a nondegenerate solution of LCP{e, M), we can use it to compute, in 
strongly polynomial time either a solution for LCP{q, M) or a secondary ray for ELCP{e,q, M); 
thus showing that the constructed bimatrix game indeed provides a 'solution' for LLCP{e,q,M). 

In Section [9] we extend the results of the previous section to accommodate the reduction of 
any LLCP{d,q,M) for which d > and LCP{d, M) is nondegenerate; thereby achieving our goal 
of showing that any LCP{q, M) which is Lemke(d) WAD-venfied (satisfying our nondegeneracy 
assumption as stated above) can be reduced to a symmetric 2-NASH. We note that for any given 
M and q, the reduction is workable for all positive covering vectors with the exception of a finite 
number of sets of measure 0. 

The constructed reduction is particulary useful since it provides a bijection between the reducible 
LCPs and their corresponding 2-NASII problems. In particular, the simplicity of the reduction and 
its bijection property allows for the practical use of the results of the extensive research on 'non 
Lemke type' 2-NASH algorithms for solving (or enumerating the solutions of) reducible LCPs . In 
addition, these reductions can be applied to investigate properties of solutions of reducible LCP 
via known properties of the associated 2-NASH problems. We discuss these subjects together with 

''a ray of ELCP{d, q, M) is an unbounded edge of ELCP{d, q, M) which includes its endpoint (a vertex of 
ELCP{q, M)) together with a direction vector. 
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additional concluding remarks in Section [TOl 

Throughout the paper we denote by e vectors all of whose entries are 1. Given a matrix A, we 
denote by j4j. the i-th row of A, by A,j the j-th column of A, and by Aij the ij-th entry of A. We 
denote by 7?-"^^"', T^-^J^^", and iV^^^ the space of m x n real matrices, the space of nonnegative 
m X n real matrices, and the space of positive mx n real matrices, respectively. Whenever n = 1 we 
abbreviate 7^"^^" to TZ^, and whenever m = n = 1 we abbreviate 7^"^^"- to TZ. 



2 LCP and Lemke's algorithm 



Given M G 72."^^"^, q G 7?."^, the linear complementarity problem, LCP{q, M), is defined as 

find z G TV^ such that 

q + Mz>0, z>0, (la) 
zT(g + Mz) = 0. (lb) 

Note that ((Ia])-([Tb]) imply 

Zi{qi + Mi,z) = 0, i = l,...,m. (Ic) 

We denote by FEA{q, M) the set of all z satisfying (fTa|) . and by SOL{q,M) the set of all z 
satisfying (fTa|) and ()lbp . 



In this section we present the generic Lemke algorithm (the so-called Scheme I - see iCPS92| . 
4.4.5). Given LCP{q, M) we define its extended version, with a covering vector d > as 

ELCP{d, q, M) = {zq G 7?.+ , z G IV]^ \ q + dzQ + Mz>Q and z'^ [q + dzo + Mz) = 0}. 

Note that ELCP{d, q, M) is composed of a polyhedral set intersected with one nonlinear com- 
plementarity constraint. Throughout the paper whenever we refer to vertices, edges and rays 
of ELCP{e,q,M) we mean the vertices, edges and rays of the polyhedral set associated with 
ELCP{e,q,M). We assume that ELCP{d,q,M) is nondegenerate, that is that the polyhedral 
set associated with it is nondegenerat^. Let {zq, z) G ELCP{d, q, M),w = q + dzo + Mz, and let k 
be the the number of positive entries in {zo,z,w). By the nondegeneracy assumption, k is equal to 
either m (in which case {zq, z) is a vertex of ELCP{d, q, M)), or m -|- 1 (in which case it is a point 
on an edge of ELCP{d,q, M)). If a vertex of ELCP{d,q, M) is contained in an edge, we say that 
the vertex is an endpoint of the edge. If an edge of ELCP{d, q, M) is unbounded then it corresponds 
to a ray of ELCP{d, q, M), which can be presented as 
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where 



(zo, z) is a vertex of ELCP{d, q, M), 
u G SOL{duQ, M), uq G {0, 1} and {uo,u) ^ 0, 
z'^{duo + Mu) = and {q + dzo + Mz) = 0. 



*There is no loss of generality in this assumption since if it is not satisfied, we perturb q by applying standard linear 
programming techniques. 
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Consider the ray of ELCP{d,q, M) with z = 0, zo = — mini<i<m uq = 1, u = 0. We call 
this ray the primary ray, and its corresponding endpoint vertex the initial vertex. Any other ray of 
ELCP{d, q, M), can be characterized as 

{zo, z) is a vertex of ELCP{d, q, M), (4a) 
u G SOL{duo, M) \ {0}, uo G {0, 1} and e^u = 1 whenever uq = 0, (4b) 
z^{duo + Mu) = {) and {q + dzQ + Mz) = (4c) 

We call such a ray, a secondary ray. We denote the set of all secondary rays of ELCP{d, q, M) by 

SR{d, q, M) = {(zo, z, uq,u) satisfying (gal-die])}. 

Note that a secondary ray {zq,z,uq,u) has two components, a vertex {zq,z) of ELCP{d,q, M) and 
what we call a secondary direction {uo,u) as defined in ()4bp . Specifically, we denote the set of all 
secondary directions of ELCP{d, q, M) as 

SD{d,M) = {(no,n) satisfying (libl)}. 

We distinguish between two types of secondary directions (and rays), according to whether uq = 0, 
which we call a type secondary direction, or no = 1, which we call a type 1 secondary direction. 
Specifically, for k = 0,1, we denote the set of all type k secondary directions of ELCP{d,q, M) by 

SDk{d,M) = {{uo,u) € SD{d,q,M) \ uq = k}. 

Similarly, for A; = 0, 1, we we denote the set of all type k secondary rays of ELCP{d, q, M) by 

SRk{d, q, M) = {{zq, z, uq, u) € SR{d, q, M) \ uq = k}. 

Starting with the initial vertex of ELCP{d,q, M), the generic Lemke((i) algorithm traces a uniqu^ 
finite path of adjacent vertices of ELCP{d, q, M), terminating with either a solution to LCP{q, M) or 
with a secondary ray of ELCP{d, q, M). Specifically, the algorithm ends with either a vertex [zq, z) of 
ELCP{d, q, M) with zq = {so z £ SOL{q, M)), or with a secondary ray {zq, z, Uo,u) G SR{d, q, M) 
with Zq > 0. We say that Lemke(d) resolves a given LCP{q, M) if either it ends with z G SOL{q, M), 
or if the terminal secondary ray can certify that SOL{q, M) = 0. Whenever Lemke((i) resolves 
LCP{q,M) we say that LCP{q,M) is Lemke(d)-resolvable. 

Ever since the introduction of the Lemke algorithm |Lem65) . extensive research efforts focused 
on identifying classes of matrices M for which LCP{q, M) is Lemke((i)-resolvable for all q. In the 
following we discuss two major groups of matrices containing almost all known classes of matrices 
M for which LCP{d, q, M) is Lemke((i)-resolvable for all q. 

The first group is based on the idea that if SR{d, q, M) = 0, then Lemke((i) outputs z G 
SOL{q,M). Specifically, we consider the class of d-regular matrices (see |CPS92j . 3.9.20) as de- 
fined below. 

Definition Given M G Tf^X"^ and d G we say that M is d-regular if SOL{dT,M) = {0} 

for all T G 72.-I-. We denote the class of d-regular matrices by R((i). 

^The uniqueness is due to the assumption that ELCP{d, q, M) is nondegenerate. 
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It follows that if M € R(d), then for all ELCP{d, q, M) has no secondary directions and thus 
no secondary rays. That is, SR{d, q, M) = for all q. Recalling that Lemke((i) terminates in a finite 
number of steps with either z G SOL{q, M) or with {zq, z,uq,u) G SR{d,q, M), we conclude that 
whenever M G R((i), LCP{q, M) is Lemke(d)-resolvable for all q. 

Remarl<|2| It is well known that M belongs to the strictly semimonotone matrix class (which is 
denoted by E) if and only if SOL{q,M) = {0} for ah g > (see |CPS92) .3.9.11). Thus, it follows 
that for all d > 0, E C Ti{d) . In addition, E properly contains the strictly copositive matrix class 
(C), and the class of all matrices whose principle minors are positive (P). Thus, we observe that 
LCP{q, M) with M in E, C or P is Lemke(d)-resolvable for all d > and ah q. 

The second group includes classes of matrices for which SR{d, q, M) ^ implies that SOL{q, M) = 
0. Specifically, most matrix classes with this property that have been identified in the LCP literature 
share the following property: 



We denote by USR((i) (for useful secondary rays) the class of all matrices M for which (I5ap -( l5bl) are 
satisfied for all g > . 

The class of matrices that satisfy (|5ap is defined below. 



Definition Given M G T^.^^^xm ^ ^ '^++' we say that M G q^rnxm -g d-semiregular if 
SOL{d,M) = {0}. We denote the class of d-semiregular matrices by Ro((i). 



1. While the term 'd-semiregular' is introduced here for the first time, the class itself has been 
introduced in |Gar73j under the name E*(d). 

2. If M G USR((i) then the existence of a secondary ray for ELCP{d, q, M) implies that SOL{q, M) = 
0. Hence, any LCP{q, M) with M G USR(d) is Lemke(d)-resolvable for ah q. 

3. R((i) C USR(d). 

4. It is well known that M belongs to the semimonotone matrix class (which is denoted by Eq) if 
and only if SOL{q, M) = {0} for ah g > (see (CPS92], 3.9.3). Thus, it follows that Eq C Ro(d) 
for all d > 0. In addition, Eq properly includes the copositive matrix class (Co), and the class 
of all matrices whose principle minors are nonnegative (Pq). 

5. There are two well known classes of matrices, L and Qo H Pq, which are known to be in 
USR(d) for all d > 0. In particular, major matrix classes, including Column Sufficient (CSU), 
Row Sufficient (RSU), and Sufficient (SU), are subsets of Pq n Qo, while Copositive Plus 
(C^), and Copositive Star (Cq) are subsets of L. Hence, LCP{q, M) where M belongs to 
any of these classes of matrices is Lemke(d)-resolvable. For a discussion of these and other 
Lemke(d)-resolvable classes see |CPS92] and |Mur88] . Figure 1 at the end of Section [5] depicts 
the relationship among these classes. 



SDi{d,M) = 

SRQ{d,q,M) ^% ^ FEA{q,M) =%. 



(5a) 
(5b) 



Remarks 



^The definitions of all the matrix classes which are mentioned in this paper can be found in [CPS92] . 



7 



3 Bimatrix Games 

Let A,B 'j^mxn ^j^g j,Qg|. niatrices of the row and column players of a bimatrix game. A Nash 
equilibrium of this game is a pair of vectors x G TZ^ , y € TZ.^ (representing mixed strategies for the 
row and column players respectively), satisfying 

Ay > e{x^Ay), B^x > e{x^By), e^x = e^y = 1, x > 0, y > 0. 

To simplify the presentation we restrict our attention to symmetric bimatrix games where A = B"^. In 
particular, it has been shown in the seminal paper [Nas51] that every symmetric bimatrix game has a 
symmetric Nash equilibrium (that is, a Nash equilibrium where x = y). In addition, it is well known 
that the Nash equilibria for any bimatrix game with cost matrices A, B (which can be assumed, 
without loss of generality, to be positive) can be easily extracted from the symmetric equilibria of 

the symmetric bimatrix game with cost matrix 

Given C G T?."'^"', we denote by SG{C) the symmetric bimatrix game where the row and column 
players' cost matrix is C. We say that x G is a symmetric Nash equilibrium of SG{C) if 

Cx > e(xTCx), (6a) 
X > 0, (6b) 
e'^x = 1. (6c) 

Note that since x'^Cx = X^^^i Xi(Cj.x), (|6al) - (l6bj) imply 

Xi{Ci,x — x^Cx) = 0, i = l,...,n. (9d) 

We denote by SNE{C) the set of symmetric Nash equilibria of SG{C). We refer to the problem of 
finding a symmetric Nash equilibrium for SG{C) as solving SG{G). 

There are several ways of formulating the problem of finding a Nash equilibrium of a bimatrix 
game as a linear complementarity problem ( |CD68| . |Eav71| . [MZ91] . |Sav0 6] ) . Here we adopt the 
reduction in jSav06j , where the problem of computing a symmetric Nash equilibrium of a symmetric 
bimatrix game is presented as a linear complementarity problem. In particular, let G be the cost 
matrix of a symmetric bimatrix game. Without loss of generality we can assume (by adding a 
sufficiently large constant to all the entries of G) that C > 0. Solving SG{G) with C > can be 
reduced to LGP{—e,G) as described in |Sav06j . and presented in the following theorem. 

Theorem 1 Suppose C > 0. 

(i) Let z G SOL{-e,G). Then, G SNE{G). 

(ii) Let X G SNE{G). Then, x^ G SOL{-e,G). 



Proof. (i) and (ii) can be easily verified by substitution. 



□ 
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4 The Complexity Class PPAD 

The class WAV [Polynomial-time Parity Argument Directed), which was introduced in the seminal 
paper |Pap94| , is a class of problems which can be presented as follows. 

Definition Given a directed graph with every node having in-degree and out-degree at most one 
described by a polynomial-time computable function f{v) that outputs the predecessor and successor 
of a node v, and a node s (which we call the initial source node) with a successor but no predecessor, 
find a node t ^ s which is either a sink (a node with no successor) or a source (a node with no 
predecessor), but not both. We call such a graph the WAT> graph associated with the problem. 

Many important problems, such as the Brouwer fixed-point problem, the search versions of Smith's 
theorem, the Borsuk-Ulam theorem and, as previously discussed, Nash equilibrium of bimatrix game, 
belong to this class |Pap94| . Interestingly, the problems in WAT) are generally believed not to be 
A/'P-hard since it has been shown in |MP91j that if there exists a WAV problem which is MV-hsid 
then NV = CoMV . What makes the study of this class attractive is that it has been shown that 
several problems within the class (such as the Brouwer fixed-point problem) are PP^P-complete 
with strong circumstantial evidence that these problems are not likely to have a polynomial time 
algorithm [HPV89j . 

The WAD complexity class seems to be a natural framework for analyzing the computational 
complexity of Lemke(d)-resolvable LCP{q, M), as the underlying graph of Lemke{d) whose nodes 
correspond to the vertices and edges of ELCP{d, q, M) has a structure reminiscent of a WAD graph. 
In particular, given ELCP{d, q, M), we define its associated graph (which we call the Lemke(d) graph 
associated with LCP{q, M)), as the directed graph G{d, q, M) whose nodes correspond to the vertices 
and edges (including rays) of ELCP{d, q, M). There is an arc {u, v) of G{d, q, M) if and only if either 
u corresponds to a vertex of ELCP{d,q, M), v corresponds to an edge of ELCP{d,q, M) and the 
vertex corresponding to u is the tail of the edge corresponding to v; or u corresponds to an edge of 
ELCP{d,q, M), v corresponds to a vertex of ELCP{d,q, M) and the vertex corresponding to v is 
the head of the edge corresponding to v. The orientations of the edges are determined according the 
scheme presented in |Tod76| . We identify the node associated with the primary ray of ELCP{d, q, M) 
as the required special source node of a WAD graph. Given (as we assume) that ELCP{d, q, M) 
is nondegenerate, we have that every node of G{d, q, M) is incident to at most two edges, and that 
there are no isolated nodes. Thus, G{d, q, M) is a nonempty collection of simple directed paths. In 
addition, any node incident to only one other node (except for the node associated with the primary 
ray) corresponds to either a solution of LCP{q, M) or to a secondary ray of ELCP{d, q, M). Thus, if 
for a given LCP{q, M) and a covering vector d, the secondary rays of ELGP{d, q, M) can certify (in 
polynomial time in the size of LGP{q, M)) that SOL{q, M) = 0, we can conclude that LGP{q, M) 
is in WAD. Whenever this is the case, we say that LGP{q, M) is Lemke(d) W AD -verified. 

Indeed, in |Pap94] , one of the first examples of a WAD problem is an LCP{q, M) where M G P. 
While it is customary in the literature of linear complementarity to discuss methods for solving 
LCP{q, M) under the assumption that M possesses some special properties, it creates difficulties 
from an algorithmic complexity point of view, as verifying these properties may be by itself a hard 
problem (e.g. identifying a P matrix is CoMV complete |Cox73| ) . Thus, in |Pap94| , the problem at 
hand (which is called P — LCP) is defined as follows. Given g,M, either find z G SOL{q, M), or 
provide a certificate (with size polynomial in the size of the problem) for M P. Motivated by the 
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discussion in |Pap94| we consider the following generic problem. 



Y — LCP{q, M) : Given M G 7?,"^^"^, q £ TV^ and a matrix class Y, find one of the following: 
(1) z € SOL{q,M), (2) a certificate that SOL{q,M) = 0, (3) a certificate that M Y. 

We say that Y — LCP{q,M) is Lemke(d)-PPAD-verified if any {zq, z,uo,u) e SR{d,q, M) leads 
(in polynomial time in the size of ELCP{d,q,M)) to (2) or (3) above. 

Remark Following the discussion in section[2l we have that USR((7, M) — LCP{q, M) is Lemke(d) 
7-''P.4I5-verified and for all q and d > 0. In particular, considering the remarks at the end of Section 
[21 we can conclude that L U (Pq H Qo) — LCP{q, M) is Lemke(d) V'PA'D-veniied and for all q and 
d>0. 



5 The generic Lemke((i) linear complementarity problem 

As stated in the introduction, it has been established that the problem of finding a Nash equilibrium 
for a bimatrix game is V'PA'D-complete. Moreover, since solving any bimatrix game is polynomially 
reducible to solving a symmetric bimatrix game, we have that the problem of finding a symmetric 
Nash equilibrium for a symmetric bimatrix game, as presented in Section [3l is also T^'P^X'-complete. 
In particular, it is shown there that this problem can be represented as an LCP{—e,M) where 
M > 0. Since M > implies that M G C (the class of all matrices for which ^ x € 7?-^ implies 
that x^Mx > 0), and considering the remark at the end of the previous section, we conclude that 
C — LCP is 'P'P^P-complete as well. In Figure 1, we display the relationship among the classes of 
matrices discussed in previous sections. An arrow from a class X to a class Y indicates that X C Y. 
So for any class Y reachable by a directed path from class C in Figure 1 we have that if Y — LCP 
is in WAT) then it is "PP^P-complete. Note that the class USR((i) (for any d > 0) contains all 
the classes of matrices Y identified in the previous section as a classes for which Y — LCP{q, M) is 
Lemke((i) VPA'D-veiiiied and for all q and d > 0. 



\ 



s 




/ 



\ 



/ 
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Next we show how Lemke((i) VPA'D-verified hnear complementarity problems can be reduced 
simply and directly to a symmetric bimatrix game. To achieve this goal, we shall consider the 
following generic problem which we call Lemke(d)-LCP(q,M) and denote by LLCP{d,q, M). 

LLCP{d,q,M) : Given M G 7^"^X"^, q G 7?."^, and d G n^^^, find either z G SOL{q,M) or 
{zo,z,uo,u) G SR{d,q,M). 

Obviously LLCP{d,q, M) is verified to be in WAV by G{d,q,M). In the following we shall 
show how it is possible to reduce any LLCP{d,q,M), whenever LCP{d, M) is nondegenerate, to a 
symmetric bimatrix game. In addition, given d > 0, we shall show, by standard LP perturbation 
techniques, that our reduction works for all d in a sufficiently small neighborhood of d, and that the 
reduction works for all d G '^++, except for a finite number of subsets of measure zero. We shall 
present our reduction in several steps, where we address the reduction of instances of the following 
problem: 

Y - LLCP{d,q,M) : Given M G U'^x^, q g U"^, d G 72-!^^ and a matrix class Y, find one of 
the following. 

(1) z G SOL{q,M), (2) {zo,z,uo,u) G SR{d,q,M), (3) A certificate that M Y. 

We start by presenting in Section[6]a simple reduction which is applicable to R(e)— LLCP(e, q, M). 
The certificate that we obtain from the bimatrix game whenever M R-(e) is of the form 7^ « G 
SOL{euo, M). So we get either a solution for LCP{q, M) or a secondary direction for ELCP{e, q, M). 
Note that at this stage we address only e as a covering vector, and that we do not reach yet our 
goal as the reduction may produce secondary directions rather than secondary rays. In Section [TJ 
we extend the previous reduction to handle Ro(e) — LLCP{e,q, M). The certificate that we obtain 
from the bimatrix game whenever M ^ Ro(e) is of the form ^ -u G SOL{e, M). So we get either 
a solution for LCP{q, M), a secondary ray for ELCP{e,q, M), or a type 1 secondary direction for 
ELCP{e,q, M). Next, in Section [HI we present a complete reduction of LLCP{e,q, M) under the 
assumption that LCP{e, M) is nondegenerate if SOL{e, M) ^ {0}. Finally, in Section [9l we extend 
the reductions in the previous sections to a general covering vector d > for which LCP{d, M) is 
nondegenerate if SOL{d,M) / {0}. 



6 Reducing R(e) — LCP{q, M) to a symmetric bimatrix game 

In this section we present a simple direct reduction of R(e) — LCP(e, q, M). In particular, given g, M, 
we construct a symmetric bimatrix game whose symmetric Nash equilibrium points correspond one- 
to-one to either z G SOL{q, M) or a certificate for M ^ R(e) in the form of 7^ u G SOL(euQ, M) 
where uq G {0,1}, so {uq,u) G SD{e,M). 

Given LCP{q, A4) with M G TZ^^^,q G 7?."^ and a covering vector e, we set n = m + 1 and a 
symmetric bimatrix game whose cost matrix C{q,M) is 

(7) 



C{q,M) 



M 




q + e 
1 
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Given C{q, M) as above, we denote any symmetric equilibrium point x G SNE{C) as x 
where y G TZ^ and t G 7?.. Given SNE{C{q, M)), we partition it to 



HI 



and 



SNE+{C{q,M)) 



SNEo{Ciq,M)) 



y 
t 



y 
t 



G SNE{C{q,M)) \ t>0 



e SNE{C{q,M)) \ t = 0}, 



In the next theorem we estabhsh a one-to-one correspondence between the symmetric Nash equihbria 
of SG{C{q, M)) which use with positive probabihty the last column of C{q,M), and the set of 
solutions to LCP{q, M). We follow this with a theorem that establishes a one-to-one correspondence 
between the symmetric Nash equilibria of G{C{q, M)) which are not using the last column of C{q, M) 
and the secondary directions of ELCP{e, q, M). 

Theorem 2 



(i) Given x 



y 

t 



G SNE+{G{q, M)), let z = y\. Then, z G SOL(ci, M). 



(ii) Given z G SOL{q,M), let t 



eTz- 



:j, y = zt. Then, x 



y_ 

t 



G SNE+{G{q,M)). 



Proof. Throughout the proof we denote G{q, M) by C. 

(i) Since t > 0, then, by ([Ml) (with i = n), t = x'^Cx. Thus, by My + {q + e)t > 

et, y > 0. Dividing by t, we get Mz + g > 0, z > 0. In addition, by (I9dp (for i = 1, . . . , m), 
yi{Mi,y + qii +i —t) = 0, so dividing by substituting for z, and summing over m, we get 
= TT=l UM.iZ + qi) = z^Mz + q). 

(ii) By (fTa]) . and setting t = ^yi+T) V — we have 

' M q + e 
1 

and obviously {e'^z)t + t = (e^z -|- l)t = 1. In addition, 

x'^Gx = y'^{My + qt + et) + t^ = t^{z^{Mz + q) + e) + l) 

Thus, since by ([Tb]), z'^{Mz + g) = 0, xTCx = ?(zTe + I) = t. Hence, x = | ^ satisfies 
(I6a])~(l6c]), and since t > 0, we have x G SNE+{C{q, M)). □ 




Theorem 3 
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(i) // 



y 





e SNEo{C{q,M)), then y € 50L(er, M)/{0} for some f > 0. 



(ii) Let u €z SOL{eT, M)/{0} for some T > 0. Then, setting y = u-^, x 



y 





G SNEo{C{q,M)). 



Proof. Throughout the proof we denote C{q, M) by C. 

(i) By dSID-debl), My > e{y^My), > y^My, y > 0, and e^y = 1. Setting f = -y^My, we get 

My + er > 0, 7^ y ^ 0- Moreover, since e^y = 1, we have (by ()9dp ) that y'^{My + er) = 0, 
concluding that y G SOL{ef, M)/{0} for some r > 0. 

(ii) Noticing that u and by p^ ~ pb|) . 

My >ei-f), yJ (My + ei-f)) = 0, y>0. 



Thus, 



y 




satisfy (|6ap - ()6bp . Noticing that e^u = 1, completes the proof. □ 



Given LCP{q, M), and combining Theorems [2] and [3l we can construct a symmetric bimatrix game 
where any symmetric Nash equilibrium point corresponds to either a solution for LCP{q, M), or a 
secondary direction for ELCP{e,q, M). Specifically, given q ^ and M, consider the symmetric 

y 



bimatrix game whose cost matrix is C{q,M). Let 
cost. We then conclude that: 



t 



G SNE{C{q, M)), and let r be its expected 



1 t > 0. By Theorem [l(i), iy G SOL{q,M), 

2 t = 0. By Theorem [1(1), y G SOL{ef, M) \ {0} for some f > 0, 

1.1 f = 0. Then, y G 5OL(0,M) with e'^y = 1, so y G SDo{d,q,M), 

1.2 T > 0, Then, yi G SOL{e,M), so y4 g SDi{e,q,M). 

So the symmetric bimatrix game SG{C) generates either a solution for LCP{q, M) or a secondary 
direction for ELCP{e,q,M). 

Remarks 



1. Note that R(e) - LLCP{e,q,M) is also reduced to SG{C{q,M)). 

2. The class of all matrices M for which LCP{q, M) is guaranteed to have a solution for all q is 
called Q. The largest known class Y which is contained in Q and for which it is known that 
Y — LCP{q, M) is Lemke(e)-resolvable, is R(e). 

3. Since P, C C E C R(e), the reduction is applicable to Y - LCP{q, M), where 1^ is P, C or E. 
Note that LCP{q, M) has a unique solution for all q if and only if M G P, and that LCP{q, M) 
has a unique solution for all (7 > if and only if M G E. 
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7 Reducing Ro(e) — LLCP{e, q, M) to a symmetric bimatrix game 

In this section we consider Ro(e) — LLCP{e, q, M) which brings us closer to achieving our goal of 
reducing any LLCP{e, q, M) to a symmetric bimatrix game. In particular, given q, M, we construct a 
symmetric bimatrix game whose symmetric Nash equilibrium points correspond one-to-one to either 
the solutions of LCP{q, M), the secondary rays of ELCP{e,q,M), or certificates for M Ro(e) in 
the form of ^ -u G SOL{e, M) which correspond to a type 1 secondary directions of ELCP{e, q, M). 

For that purpose we introduce the augmented problem LCP{q, M) associated with LCP{q, M), 
where 



M = 



' 1 -eT ■ 






e M 


, Q = 


. 1 . 



and (3 > e^ z for any vertex (not necessarily feasible) (zo)-^) of ELCP{e,q, M). 
Remarks 



1. It is a standard result in LP theory that if the entries in q,M are rational then f3 is of size 
polynomial in the size of LCP{q, M), and that (3 can be computed in time polynomial in the 
size of LCP{q,M). 



2. Augmented LCP systems where Mn is equal to (see |CPS92j ) or —1 [Tod73] are used in the 
LCP literature to eliminate secondary rays. Such augmentations do not work in our case since 
the reduction of LCP{q, M) to a symmetric bimatrix game would yield a pure Nash equilibrium 
(using with probability 1 the strategy corresponding to the first column of M) which yields 
no information about the solution (or lack thereof) of the original LCP{q, A-I). To avoid this 
possibility, we need Mn > 0, hence the choice of 1. 

In the following theorem we establish the relationship between LCP{q, M) and LCP(q, M). 
Theorem 4 



(i) 

(ii) If 



E SOL{q, M) if and only if z e SOL{q, M). 

€ SOL{q,M) where zq > 0, then there exists {zq,z,Uo,u) G SR{e,q, AI) and A > 



z 

such that 







zq 


+ A 


no 


z 




z 




u 



Then there exists A such that 


Zo 


+ A 


no 


G 




z 




n 





G SOL{q,M). 



Proof. 



(i) The 'only if direction is obviously true. The 'if direction is true because of the nondegeneracy 
assumption (so 5 is a a vertex of LCP{q, M)) and by the definition of /3. 
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(ii) Let 



z 



e SOL{q, M) where zq > 0. Then, 
q + ezo + Mz >0, z>0, zq > 0, z'^{q + ezo + Mz) = 0, 5o(/3 + zq - e^z) = 0, 



which imphes that {zq,z) € ELCP{e,q, M) and (since zqj > 0) e^z = /3 + ig- However, by the 
definition of /3, (zo,z) must be a point on a secondary ray of ELCP{e,q, M). That is, there 



satisfies all the 



exists A > and {zq, z,uq,u) G SR{e,q, M) such that 
(iii) By the definition of SR{e, q, M), we have that for all A > 0, 







zo 


+ A 


Uo 


z 




z 




u 



zo 


+ A 


no 


z 




u 



constraints of LCP{q, M) except possibly for the last constraint. However, since —z^ + e^z < f3, 
u / 0, no e {0, 1}, and -no + e^n > 0, setting A = yields -(zo + Ano) + eT(z + An) = /3 

which, considering that A > 0, completes the proof. □ 



Remark The extraction of either a solution or a secondary ray from a secondary direction as 
described in the proof of Theorem [3]-(ii) can be done by standard LP technique that can be executed 
in strongly polynomial time (that is, the required number of elementary calculations such as additions, 
multiplications, divisions and comparisons is bounded above by a polynomial function of m). 

Next, we show that M R.(e) implies that M Ro(e), which allows us to apply the reduction 
of the previous section to the augmented problem. 



Theorem 5 //O / 



no 
u 



G SOL{eT,M) for some r > 0, then no + r > and G SOL{e,M). 



Proof. By the premise of the theorem there exists no > 0, n > 0, no + e^n > 0, such that 



M 



no 
n 



+ 



r > 



, and 



[Uo W 



" 1 -eT ■ 




Uo 


+ 


" 1 ■ 




' " 


e M 













u 


e 



Thus, 



no — e^n + r > 0, 

Mu + e(no + r) > 0, n > 0, 

nT(Mn + e(no + f)) = 0. 



(8a) 
(8b) 
(8c) 



By (jSap we have that no + r > (as otherwise n = 0,no = contrary to the assumption). Thus, 
from dHEl) and ([HcD we have that ^^u G SOL{e, M). □ 



Combining Theorems S] and El and recalling the definition of the class Ro(e), we get that, given 
LCP{q, M), we can construct a symmetric bimatrix game where any symmetric Nash equilibrium 
point corresponds to a solution of LCP{q, M), a secondary ray of ELCP{e, q, M) or a type 1 direction 
of ELCP{e,q,M). 
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Let 



Specifically, consider the symmetric bimatrix game whose cost matrix is C{q, M) = 
G SNE{C{q, M)), and let f be its expected cost. We then conclude that: 



1 -eT -/3 + 1 
e M q + e 
1 



z 

t 



1 t > 0. By Theorem m 



z 



G SOL{q,M). 



1.1 zq = 0. Then (by Theorem ID- (i)), z~ G SOL{q, M) (a solution to the original problem). 

1.2 Zq > 0. Then (by Theorem |3]-(ii) and the remark following its proof), we can obtain (in 
strongly polynomial time) (zq, z, uo,u) G SR{e, q, M) (a secondary ray of ELCP{e, q, M)). 



2 t = 0. By Theorem [3 / 



no 
u 



G SOL{eT, M) where r > 0. Thus, by Theorem [5l 



^ ^li? e SOL{e,M) (so 
ELCP{e,q,M)). 



G SDi{e,q, M) (a type 1 secondary direction of 



8 Handling nondegenerate type 1 secondary directions 



In this section we show that if 7^ n G SOL{e, M) is nondegenerate, then we can compute, in strongly 
polynomial time, (zq, z), a vertex of ELCP{q, M) such that if zq > then (zq, z, 1, u) G SRi{e, q, M) 
(that is, a type 1 secondary ray). 

Let M be a non-zero, nondegenerate solution for LCP{e, M). Setting v = Mu + e, let a = {i\ui> 
0} and a = {i \ Vi > Note that by the nondegeneracy assumption a U q = {1, . . . , m} and Maa 
is nonsingular. Now, we set 



Za 




' M„„ ■ 


-1 


qa 




Za 




' " 


Wa 












Wa 








which results in 



w 



+ Mz, z'^w = 0, u^{q + Mz), z^{Mu). 



Thus, if z,i/) > then z G SOL{q,M). Otherwise, since Uq > and Va > 0, for sufficiently large 
A > we have that 



Za 
Wa 



+ A 



> 








Letting A be the smallest A satisfying the inequality above, and setting zq = X, z = z + Xu we get 
that {zq,z,1,u) G SRi{e,q, M). Note that constructing (zo,z) whenever a nondegenerate 
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9 Extensions 

In this section we show how to extend our reductions whenever a general positive covering vector is 
used rather than e. The key to the results in this section is the following proposition. 



Proposition 6 Given d G ^ diagonal matrix with Da = di. Then, (^0,2:) G 

ELCP{d,q,M) if and only if{zo,Dz) £ ELCP{e, D'^q, D'^ MD-^). 

Proof. The proposition is easily verified by observing that 

q + dzo + Mz>0, z>0 if and only if D-^{q + dzQ + MD'^Dz) >0, Dz> 0, 

and zT(g + dzo + Mz) = if and only if z^DD-^{q + dzo + MD-^Dz) = 0. □ 

Corollary 7 Given d € TZ^^, let D be the diagonal matrix with Da = di. Then, 

(i) z G LCP{q,M) if and only if Dz G LCP{D-^q, D-HlD~^). 

(ii) M G R(d) if and only if D-^MD-^ G R(e). 
(ill) M G Ro((i) if and only ifD-'^MD-'^ G Ro(e). 

Proof. 

(i) Results from Proposition [6] by considering zq = 0. 

(ii) -(ii) Result directly from (i) and the definitions of R((i) and Ro((i). 

By Proposition [6] and Corollary [7] it can be readily verified that the results of Section [6] can be ex- 
tended to K{d)-LCP{q, M) by considering B^{e)-LCP{D-^q, D-^MD-^). Similarly, the results of 
SectionElcan be extended to Yio{d)-LLCP{d, q, M) by considering Ro(e)-LLCP(e, D'^q, D-'^MD'^). 
Finally, given that LCP{d, M) is nondegenerate we can (by PropositionEl Corollary[71 and the results 
of sections [7] and [8|) reduce any LLCP{d, q, M) where LCP{d, M) is nondegenerate to a bimatrix 
game whose cost matrix is 



C 



1 -PY{T=idi + l 

e D-^MD~^ D-^q + e 
1 



Let 



G SNEiC). Based on Proposition [6] and Corollary [7] we replace z with D z and proceed 



z 
i 

to recover a 'solution' to LLCP{d.q.M) by following the steps prescribed in sections [7] and [HI 



17 



10 Concluding remarks 

1. The main result of this paper is that for almost any given d € 'R-^_i^,q € TZ"^, M E 72."*^"^, it 
is possible to effectively set up a symmetric bimatrix game whose Nash equilibria correspond 
one-to-one to all the endpoints (excluding the one corresponding to the primary ray) of the 
directed graph associated with LCP{q, M) and the Lemke method with a covering vector d. 
The only condition is that any Nash equilibrium corresponding to a type 1 secondary direction 
(that is a solution for LCP{d, M)), has to be nondegenerate. Note that if this is not the case, 
we can perturb (by standard LP techniques) d to a d which, when used as the covering vector, 
will guarantee that the reduction will work. This observation also means that for any given 
M, q the reduction is workable for all covering vectors d E _(- 1 with the exception of a finite 
number of sets of measure 0. 

2. As a consequence of the main result as specified above, the reduction will resolve any LCP{q, M) 
which is Lemke((i) "PP^P- verified. Note that all the major matrix classes of M which are 
known to be Lemke(d) resolvable for all q are actually known to be Lemke(o?) PT'^P-verified 
(and typically for all d G 7?.^_|_). These classes (which are all subsets of USR((i)) and their 
relationships are depicted in Figure 1. 

3. The direct reductions which are presented in this paper highlight the importance of the problem 
of 2-NASH within mathematical programming. In a sense, we show that as any LP can be 
directly reduced to zero-sum game (see |Dan51j and [Adll2|), it is analogously possible to 
directly reduce any LCP{q, M) which is Lemke(d) P'P^X'-verified to a 2-NASH problem, thus 
showing that many of the results regarding such problems are relevant to LCP theory. 

4. The reductions in sections El El andEjare simple and easy to execute. Thus, any algorithm that 
is applicable to bimatrix games can be directly used to solve instances of RSU((i) — LCP{q, M) 
and LLCP{d,q, M). Also, considering that the proposed reductions are bijections, algorithms 
with a variety of goals, such as enumerating all, or specific subsets of Nash equilibrium points 
can be applied for similar goals regarding the solutions of linear complementarity problems for 
which our reductions are applicable. It should be noted that there is a vast literature covering 
the subjects of computing and enumerating Nash equilibria of bimatrix games (see e.g. the 
surveys in jvSt02] . |vSt07] and the papers introduced in [vStlO| ). 

5. Over the years several refinements of Nash equilibrium have been introduced. In particu- 
lar, some results regarding the existence and computation of these refinements have been 
established. In |MT98j some of these refinements are generalized to LCPs. The reduction 
of LCP{q, M) which are Lemke{d) VPA'D-vehiied (e.g. where M G USR((i)) to symmetric 
NASH-2, provides us with a tool to investigate analogous questions with respect to the general- 
ized refinements to such LCPs. For example, in a forthcoming paper, we demonstrate such an 
analysis by proving that any LCP{q, M) with M G R(e) has a proper solution. As a corollary 
of this analysis we prove that the (unique) solution of LCP{q, M) where M G P is proper and 
thus settle a conjecture posed in |MT98| (where it is proved for 2x2 matrices). 

6. The simple reductions proposed in this paper allow us, whenever applicable, to potentially 
gain additional insight into the nature of models represented by these LCPs. This seems to be 
especially useful for economic models such as market equilibrium. 
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